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Abstract. We study spectral properties of Schrodinger operators on 
R"^. The electromagnetic potential is assumed to be determined locally 
by a colouring of the lattice points in Z"^, with the property that fre- 
quencies of finite patterns are well defined. We prove that the integrated 
density of states (spectral distribution function) is approximated by its 
finite volume analogues, i.e. the normalised eigenvalue counting func- 
tions. The convergence holds in the space L^{I) where I is any finite 
energy interval and 1 < p < oo is arbitrary. 



1. Introduction 

Spectral properties play a key role in the analysis of selfadjoint operators. 
This is in particular the case for Hamiltonians describing the time evolution 
of quantum mechanical systems. In the context of mathematical physics one 
often studies the integrated density of states, in the following abbreviated 
by IDS. It is very natural to think of the IDS as the normalized eigenvalue 
counting function of the restriction of the Hamiltonian to a large but finite 
volume system. This leads to the question in what sense and how well one 
can approximate the IDS of the full Hamilton operator by the spectral dis- 
tribution functions of appropriately chosen finite-volume analogues. This 
question has been pursued for various types of selfadjoint operators, resp. 
Hamiltonians, in particular in the mathematical physics and geometry lit- 
erature. Let us mention the seminal papers |Pas7H IShu79| and the recent 
reviews |KM071 IVesOS] . There one can find also an overview of the literature 
up to '07. 

It turns out that in the discrete and in the one-dimensional setting one 
can control the convergence of finite volume approximants to the IDS very 
well. Let us state this more precisely: 

• For difference operators (with finite range) on combinatorial graphs 
the eigenvalue counting functions are bounded. This leads to uniform 
convergence (in supremum norm) for the IDS |LMV08| ILV09] . 

• For Schrodinger operators on metric graphs (so called quantum 
graphs) with constant edge lengths one can achieve uniform con- 
vergence for the IDS as well |GLV07j . Here one has to assume that 
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the randomness satisfies a finite local complexity condition. Note 
that for quantum graphs the eigenvalue counting functions are un- 
bounded: However, the technically relevant objects are the spectral 
shift functions, which are still bounded. 
• Metric graphs with non-constant edge lengths lead to unbounded 
shift functions; in this case, convergence holds locally uniformly, as 
well as globally uniformly with respect to a weighted supremum norm 
[GLVHS] . 

See also |GLV08j for an overview. 

For electromagnetic Schrodinger operators on M.'^ even the perturbation 
by a compactly supported potential may lead to a locally unbounded spec- 
tral shift function. Thus the shift function diverges not only at infinity but 
also on compact energy intervals. This is in particular the case for Landau- 
type Hamiltonians, see, e.g., |RW02l [HKN^06] and references therein. One 
may expect that the situation will be better for certain random perturba- 
tions of the Landau Hamiltonian. Li order to obtain continuity of the IDS of 
Landau-type Hamiltonians plus a random, ergodic potential one has to pose 
appropriate conditions on the randomness. They amount to regularity con- 
ditions on the random distribution (see |CH96l IWan97l IHLMWOll ICHK07] 
and references therein). The results of the present paper apply to highly 
"singular" distributions, though: those of Bernoulli type; at each lattice site 
in Z*^, local electric and magnetic potentials are chosen randomly from a fi- 
nite set of prototypes. For such models there are no results on the continuity 
of the IDS. Thus this property cannot help us proving uniform convergence 
of the distribution functions. Our main result is that one can still achieve a 
strong form of convergence. More precisely, we show that convergence holds 
in the space L'P(I) for any finite interval / C M and any finite p. 

In the following section we describe our model and assumptions and state 
the main theorems. Section [3] provides bounds on the spectral shift function. 
They are applied in Section 2] which establishes certain almost additivity 
properties and thus concludes the proof of the main theorem. The latter is 
applied to certain types of alloy type random Schrodinger operators in the 
final section. 

2. Model and results 

Throughout the paper we will consider electromagnetic Schrodinger op- 
erators which satisfy the following regularity 

Assumption 1. Let lA he an open set in , A: lA ^ a magnetic vector 
potential, each component of which is locally square integrable, V = V-^ — 
V-: U M a scalar electric potential, such that its positive part V^- > 
is locally integrable and its negative part V- > is in the Kato class. This 
implies that V- is relatively form bounded with respect to — A^, the Dirichlet 
Laplacian on lA, with relative bound 5 strictly smaller than one. Under these 
conditions the magnetic Schrodinger operator 



(1) 



= (-iV -Af + V 
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is well defined via the corresponding lower semi-bounded quadratic form with 
core C^{U) |Sim79j . Due to this choice of core, we say that has Dirichlet 
boundary conditions. 

Let us mention locally uniform L^-integrability conditions which are suf- 
ficient for V- to be in the Kato-class. More precisely, if V- satisfies 

ll^-ilLf .,(Rd) = sup ( / \V^{y)\Pdy) < oo 



xeM.'' ^ J \x-y\<l 

for p = 1 if d = 1 and p > d/2 if d > 2, then it belongs to the Kato-class. 

Next we want to introduce the notions of a colouring and a pattern. For 
this purpose we denote the set of all finite subsets of Z*^ by J-fin(Z'^) and an 
arbitrary finite set by A. A colouring is a map 

C:Z'^ ^A 

and a pattern is a map P : D{P) A, where D{P) G J'fin(Z'^) is called 
the domain of P. We denote the set of all patterns by V. For a fixed 
Q € J-fin(Z'^) we denote the subset of V which contains only the patterns 
with domain Q by V{Q). Given a set Q C D{P) and an element x G Z*^ we 
define a restriction of a pattern P\q and the translate of a pattern P + xhy 
the vector x in the following way 

P\q: Q ^ A,g ^ P\Q{g) = P{g), P + x: D{P) + x ^ A,y + x ^ P{y) 

Two patterns Pi , P2 are equivalent if there exists an x € Z'^ such that P2 = 
Pi + X. The equivalence class of a pattern P in "P is denoted by P. This 
induces on "P a set of equivalence classes V. For two patterns P and P' the 
number of occurrences of the pattern P in P' is denoted by 

ttp(P') := «{x G Z'^ I D{P) +XC D{P'), P%(p)+, = P + x}. 

Here, ft denotes the cardinality of a finite set. 

Next we define the notion of a van Hove sequence and of the frequency 
of a pattern along a given van Hove sequence. A sequence (C/j)jgN of finite, 
non-empty subsets of Z'^ is called a van Hove sequence if 

(2) for ah M G N : lim ^ = 0. 

Here d^U = {x € U \ dist{x , Z'^ \U) < M}U{x £ Z'^XU \ dist(x,C/) < M}. 

It is sufficient to check the relation ([2]) for M = 1, it then follows for all 
M G N, cf. for instance Lemma 2.1 in [LSVIO] . If for a pattern P and a van 
Hove sequence {Uj)ji=^ the limit 



Up := lim 



exists, we call up the frequency of P along {Uj)j,=f>} in the colouring C. 

In our setting A will be a finite collection of pairs (a, w) where a is a function 
M"^ — 7> M'^ such that all its components are in and w is a function R'* — > M 
such that its positive part is in and its negative part is in the Kato class. 
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cf. Assumption [TJ Moreover, both the support of a and of v are contained 
in 

Wo := [0, l]'^. 

Given a colouring C: Ij'^ ^ Awe denote by Ca its first and by Cy its second 
component. To each colouring C we associate an electromagnetic potential 
{Ac,Vc): M'^ ^ M'^ X R 

Ac{x) := Ca(A:)(x - k), Vc{x) := ^ C,{k){x - k) 

fcezd 

and a Schrodinger operator 

He ■■= i-iV-Acf + Vc. 

Note that for any open subset U of Mf^ the restriction of Hq to U with 
Dirichlet boundary conditions satisfies Assumption [TJ 

Now we want to define the IDS of Hq and of finite restrictions thereof. For 
this purpose we need some more notation. In order to associate finite subsets 
of Z'^ with bounded subsets of M'^, we define 

W: J-finlZ'^) ^ B{R''), Q^Wq:=\J {Wo + t) 

teQ 

where we use the natural embedding Z"^ C M'^ and denote the Borel-c-algebra 
by B. Furthermore, if U is an open set in M'^, a Schrodinger operator 

o 

defined on U and Q € J^fin(Z'^) such that Wq cU we define, by a slight but 
natural abuse of notation, 

i.e. the restriction of to the interior of Wq in the sense of quadratic 
forms as above. 

Let us denote by X(_oo,A] i^c) and X(-oo,x] {^c) spectral families of 
He and H^ , respectively. The IDS of H^ is the distribution function 

N{X,Q) :=T:[x(_oo,A](^?)], 

divided by the volume volWg = ^Q. Since Q is finite, H^ is an elliptic 
operator on a bounded domain and X(_oo,A] i^c) is trace-class. However, 
X(-oo,A] i^c) is not, which is the reason why we need the existence theorem 
below, and why Hq may display any interesting spectral features at all. 

For M G N we denote by Cm C Z"^ the cube at the origin with side length 
M- 1, i.e. 

Cm ■■={x eZ'^ ■.0<xj <M - l,j = l,...,d}. 

Now we are prepared to state our main theorem: 

Theorem 2. Let C: Z"^ — )• ^ be a colouring, (C/j)jgj^ a van Hove sequence 
such that for all patterns P G UAfeN^(^A^) frequencies up exist. Let 
I C M be a finite interval. Then there exists a function N belonging to all 
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L'P{I) with 1 < p < oo and independent of the van Hove sequence such that 
for j — )■ oo we have 



(3) 



for any any p € [1, oo). More precisely, for any Af G N the above integral is 
bounded by 



PeViCM) 



where T = sup/. The dependencies of the constants appearing in the esti- 
mate are as follows: Cp^d depends only on the dimension d and the exponent 
p, and C depends only on d and on (the Kato norm of) Vc-- 



Remark 3. • Obviously, the theorem yields a function A'^ G 

such that on each finite interval / the above holds, and we may 
extend this to (upper) semibounded intervals since our operators 
are uniformly semibounded below. 

• In particular one may choose the van Hove sequence Uj := Cj or a 
similar family of expanding cubes (if the frequencies exist) so that 
one gets an explicit i-decay for the second term in the error estimate. 
A sequence of cubes is a common choice for defining the IDS. 

• In the random setting, introduced in Section^ there is an alternative 
definition of the IDS (Pastur-Shubin formula) which coincides with 
N, as we show there. 



3. Bounds on the spectral shift function of facets 

In this section we prove certain bounds on the spectral shift function 
(SSF) which are needed in Section 21 They concern the SSFs of two elec- 
tromagnetic Schrodinger operators which differ by an (additional) Dirichlet 
boundary condition on a facet. This difference can be understood as a gen- 
eralized (positi ve) compa ctly supported potential. Hence we can use results 
established in [HKN^OG] , either verbatim or with slight modifications. 

Recall that if V- is in the Kato class, then there exists some number 5 
smaller than one such that V- is relatively A-bounded with relative bound 
6. 

We quote Lemma 5 from |HKN"'"06| : 

Lemma 4. Let be a Schrodinger operator defined on the open set lA C 
M"^ of finite volume which satisfies Assumption [7]. Denote by En the n^^ 
eigenvalue counted from below including multiplicity of . Then there 
exists a constant Ci such that 

(4) £„>M1_j5MQ^''_c, forallneK 
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The constant Ci depends on the Kato norm of V- only. Lemma U] will 
be used in the proof of the next result, which is a slight modification of 
Theorem 1 in [HKN+06] . 

We will frequently use certain subsets of d — 1 dimensional hyperplanes 
in R"^. They are unit squares in the hyperplanes. More formally we define: 

Definition 5. A set S" C M'' is called canonical facet if there exists a j € 
{1, . . . ,d} such that 

S = {(xi, . . . , Xii) \ Xj = and Xi € [0, 1] for i ^ j}. 

A set S is called a facet if there exists a canonical facet S and a vector 
X G Z'^ such that 

S = X + S := {y eR'^ \ y - X e S}. 

Let U be an open subset of M'^, S a facet as defined in Definition O and 
U := L{\S. Let Hi be a Schrodinger operator on U, satisfying Assumption[TJ 
Using quadratic forms we define H2 as the Dirichlet restriction of Hi to U. 
Then the operators e~^^, e~^^, and Vqs ■= e~^'^ —e~^^ are well defined by 
the spectral calculus. 

Theorem 6. (a) The operator Ves is compact. 

(b) Denote by the n^^ singular value ofV^Q counted from above including 
multiplicity. Then there are finite positive constants c and C2 such that the 
singular values of the operator Ves obey 

(5) fin<C2e-''^''\ 

The constant c may be chosen depending only on d, while C2 depends only 
on the Kato-class norm ofV-. 

Proof. To prove part (a), i. e. that the operator V^s is compact, it is sufficient 
to find a family Aj^ ,R > of compact operators such that the operator norm 
of the difference := Ves — Ar tends to zero as i? — >■ cx). Indeed, such an 
operator family will be constructed in the proof of the quantitative statement 
(b). The operators in the family will be even trace-class. 

The proof of (b) is almost the same as the one of Theorem 1 in [HKN+n6| . 



We use the same notation as there and explain only the step in the proof 
which is slightly different in the present setting. Let B := Br C be an 
open ball of radius R > containing the facet S in its interior, and the 
complement of S. 

Denote by Hf the Dirichlet restriction of Hi to the set UCiB and by H2 
the Dirichlet restriction of H2 to the set U B. Set 

Let Ea; and P^^ denote expectation and probability for Brownian motion, bt, 
starting at x. For any open set U CM.'^ denote by tu = inf{t > \ bt ^ U} 
the first exit time from U. We use the Feynman-Kac-Ito formula to express 
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where 5^ is a real valued stochastic process (Ito integral) corresponding 
to the magnetic vector potential A of the Schrodinger operator; this repre- 
sentation holds for more general / [Sim79j but a dense set suffices for our 
purposes. Analogous representations hold for the operators e~^^, e~^? , 
and e~^2 if one replaces the condition X{tu>i} by 

X{r^>i},X{r^,nB>i}' ^d X{r^^3>i}, respectively. 

Since the operator D can be expressed in terms of four different exponentials 
it follows that 

{Df){x) = E,. [p{h)e-'''Ai^'^e-^oy^^^)'^^f{h^) 

where 

P = X{ru>l} - X{r^>l} - X{runB>l} + X{r^^^>l}- 

A simple transformation (using X{rM>i} - X{t^>i} = X{rM>i}(l - X{t^>i}) = 
X{7i/>i}X{rjj<i} = X{tu>i}X{tsc<i} etc.) shows that 

P = X{tu>1} X{tsc<1} X{ts<1}- 

We abbreviate B := {t5c < 1} n {tb < 1}. The Holder inequality implies 
that 

|Z)/|(x)<(E4x{^>i}e-^^o^(''^)'^l)'^'(E4xM&)])'^'(E4|/(6i)P])'^'. 

At this point we make the dependence of the operator D on the radius of 
the ball B = explicit and denote it consequently by Dr. Prom this point 
on we can follow exactly the proof of Theorem 1 in [HKN"'"06] to conclude 
that 

\\Dr\\ < const exp 



32 

With the choice R = Rn = 'n}/'^'^ the desired estimate ([5]) follows. □ 

The next Lemma is an abstraction of the proof of Theorem 2 in HKN"'"06] , 
which in turn relies on |HS02] . The abstract formulation may be of use also 
in other contexts. 

Let A, B be two selfadjoint operators such that V^s = — is trace 
class. This implies that the sequence ^ = {//njnsN of singular values of V^s 
(enumerated in decreasing order including multiplicity) converges to zero. 
Let F: [0,oo) — >• [0,oo) be a convex function with -F(O) = 0. In particular, 
F is isotone because for x G [0, oo),a G [0,1] convexity gives F[ax) = 
F{ax + (1 - a)0) < aF{x) + (1 - a)F(O) = aF{x) which, since F{x) > 0, 
implies F{ax) < F{x). Set (j){n) = F{n) — F{n — 1) for n € N, n > 2. 

Lemma 7. (a) Let F he as above. Then 



(6) / F(|e(A,i3,A)|)dA<e^(</.,M)^2(i,) 

J — oo 

(h ) Let /i : M — 7- M 6e a bounded measurable function with support in (— oo, T] . 
Then 

(7) / /i(A)e(A,i?,.4)dA<e^(0,/i),2(p,)+ T G{\h{X)\) dX 
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where G denotes the Legendre transform of F, i.e. G{y) = sup{xy — F{x) \ 
X >Q} for y > 0. 

Of course tlie usefulness of the Lemma depends heavily on a priori in- 
formation about the decay of the sequence fi. However, for the choice of 
operators A = Hi and B = H2 we do know that the singular values decay 
at an almost exponential rate. 

Remark 8. Depending on how many additional properties we assume for 
the function F, we obtain correspondingly more information about its Le- 
gendre transform G. This will be discussed next. 

(1) The following properties hold under no additional assumptions on 
F: G(0) = 0, G is convex (because xy — F{x) is convex in y) and 
G{y) > for all y. 

(2) If lim = 00 then G takes on finite values only (and vice versa). 

(3) If F is twice differentiable, f := F' > and F" > on [0, 00), then 



ify</(0), 
yf-\y) - F{f-Hy)) ify >/(0). 



G{y) 

Here f~^ denotes the inverse of the function /. Consequently, 



-00 



G{\hiX)\) dX= J {\h{\)\f-\\h{X)\) - F{f~H\h{X)\))) dX. 

{\<T:\h(X)\>f{0)} 

(4) If there exist a positive constant G and an exponent p > 1 such that 
l^n < Gn~^, then one can choose the function F as F{x) = x'^'^^, 
where q is any number smaller than p — 1. Indeed, for this choice of 
F we have 4>{n) < (q + l)n'^. Thus 

^„0(n) <{q+ 1)G n-Pn" < 00. 
n n 

Note that in this case G(y) = ( j ■ 

(5) If there exist positive constants c, C and p such that /i„ < Ce"^"'' 
for all n G N, then for each value of t < c, the choice 



F{x) = r {e'y - 1) dy 
Jo 



gives a finite right hand side in Indeed, (t){n) = f^_^ (e*^'' 
^) dy < e*"'' and thus 

(<A,/i),2(i,)<C^e— ''e*-''<oo. 
n 

The Legendre transform for such a choice of F satisfies 
G{y) < yf-\y) = (^^^^.^±A^"' for all y > 0. 



LP-APPROXIMATION OF THE INTEGRATED DENSITY OF STATES 

Thus in this specific case inequahty ([7]) reads 
(8) / /i(A)e(A,S,A)dA < e^C^e~("^*)"' + 

neN 



which recovers the result of HKN"'"06 

Now we prove Lemma [71 

Proof. Since V^q is trace class but not necessarily the operator difference 
A — B, the SSF is defined via the invariance principle 

r FmX,B,A)\)dX = r F(|C(e-\e-^,e-^)I)dA. 

J —oo J —oo 

Now a change of variables gives us 

/T f'OQ 
-oo J e ^ 

To the last expression we can apply the estimate of [HS02] and bound it 
above by 

/oo 
F(|e(s,e^^,e-^)|)(is < ^/i„(V;ff)0(n). 

-"^ nGN 

This establishes claim (a). To prove (b), we note that by the very definition 
of the Legendre transform, the Young inequality 

\h-^\<Fm+G{\h\) 

holds. Integrating over A we obtain 

h{X)CiX)dX< I F{\i{X,B,A)\)dX+ j G{\h{X)\)dX. 

J — oo J 

Together with (a) this completes the proof. □ 

4. Almost additivity for the eigenvalue counting functions 

In this section we prove Theorem [2l For this aim we will apply a Ba- 
nach space- valued ergodic theorem obtained in [L MVOSj . Actually, for our 
purposes it will be convenient to quote a slightly streamlined version of this 
result from [LSVlOj . To spell it out we need to introduce the notion of a 
boundary term and the properties of almost additivity and invariance. 

Definition 9. A function b: Ffi^{'L'^) — ?> [0, oo) is called a boundary term if 
the following three properties hold: 

(i) b{Q) = b{Q + x) for ah x G and all Q £ 7'fin(Z'^), 

(ii) lim = for any van Hove sequence and 



J-+00 



(iii) there exists a constant D € (0, oo) such that 

KQ) <D^Q for ah Q G Fa^i^'^) 
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Definition 10. Let {X, || • ||) be a Banach space and F a function — )• 
X. 

(a) The function F is said to be almost- additive if there exists a boundary 
term b such that 

m m 

||F(ur=iQfc) - j;F(Qfc)|l <Y.^{Qk) 

k=l k=l 

for all m G N and all pairwise disjoint sets Qk G J-fin(^'^), k = 1, . . . ,m. 

(b) Let C : — > ^ be a colouring. The function F is said to be 
C-invariant if 

F{Q) = F{Q + x) 
whenever x € Z'' and Q € J-fin(Z'^) obey C\q + x = C\q+x- 

In this case there exists a function F defined on the (classes of) patterns 
such that F{P) = F{Q) if C\q = P. 

If F : Jfin(Z'^) — >■ X is almost additive and invariant there exists a € 
(0, oo) such that 

(9) \\F{Q)\\<K^Q for all Q e -Ffin(Z'^). 

Now we are in the position to quote the Banach space valued ergodic 
theorem from |LMV08] . see [LSVlOl §5.1] as well. 

Theorem 11. Let A be a finite set of colours, C: A a colouring 

and (Uj) a van Hove sequence along which the frequencies of all patterns 
P G UMeN'^(^Af) exist. Let F : T(JI^) X be a C-invariant and almost- 
additive function. Then the limit 

j—^oo uUn A/— >-oo uC M 

^ ^ P&V{Cm) ^ 

exists in X. Furthermore, for j,M G N the bound 



holds. 



We want to apply the ergodic theorem to the eigenvalue counting func- 
tions of Schrodinger operators, considered as elements of X := L^{I) for a 
fixed finite interval / C M and p G [1, oo). 

More precisely, we study the function 

F: J-fi„(Z^)^LP(/), Q^iV(.,/7Q) 

with N{\,H^) := Trx{-oo,X]{H'^) for A G M. Note that this notation is 
slightly different from the one used in Theorem [2j The reason is that in the 
proofs we use a more general class of operators than which was necessary to 
formulate the main result and thus need a bit more flexibility. To conclude 
Theorem [2] we need to show that F fulfils the hypotheses of Theorem [TTJ 
This is done in the following 

Lemma 12. The function F : Ffin{1j'^) — )• LP{L) is invariant and almost- 
additive. 
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Proof. Note that H is a local operator, thus depends only on A\u and 
V\ii, for any open U C M'*. The translation operator Tyf{x) = f{x — y) 
is unitary for any y € W^. Thus the spectrum of H, resp. , is invariant 
under conjugation by Ty. It follows that the function F is C-invariant. 

To prove almost-additivity, let Q = UfcLi Qk with pairwise disjoint Qk S 
F^^{l/),k = l,...,m. We need to compare F{Q) to X]fcLi-^(Qfc)- Note 
that Wq = UfcLi ^Qk iised not be pairwise disjoint because 

their boundaries can touch. There are two extreme cases: 

o o 

(1) All Wq^ are pairwise disjoint. Then Wq = UfcLi^Qfc> ^-^id conse- 
quently = 0^^, H^'^ and F{Q) = ^Li F{Qk)- 

(2) No Wq^. is disjoint from all others. Since dWq^ consists of at 
most 2d\))dQk facets (where dQk denotes the combinatorial bound- 

o o 

ary of Qk C Z'^) the sets Wq and UfcLi ^Qk differ by at most 

2c^Er=i S^Qfc facets. 
In fact, the latter case gives an upper bound for the general case (where 
the "isolated" Qk simply do not contribute), and we can drop the factor 
2 because touching facets need to be counted once only (they are counted 
twice in the sum). 

o 

So, let M < dJ2T=i ^'^Qk be the number of facets by which Wq and 

o 

UfcLi ^Qk differ and enumerate them arbitrarily as Si, ... , Sm- Set Hq := 
and 

j 

Hj := H^^ where Uj:=WQ\[j Si 

1=1 

for j = 1, . . . , M. Clearly, Hm = ®T=i and iV(-, Hm) = ET=i F{Qk)- 
This relation allows us to write the difference that we want to estimate as a 
sum of SSFs: 

m M M 

F{Q) -Y^HQk) = E i^i-^H^.i) - N{;H,)) =J2a;H,-i,Hj) 
k=i j=i j=i 

Hj^i and Hj differ exactly by a Dirichlet condition at one facet Sj so that 

Theorem [6] applies and gives the estimate /i„ < Ce"'^"^'''' for the singular 
values. Now we can apply Lemma [7| with T := sup/, A := Hj, B := Hj^i 
and F{x) := x^. Then inequality (0) reads 



(11) / \a;H,^i,H,)\P <e'^y2C2pnP-'e- 



with a constant C which is independent of j. The triangle inequality thus 
gives 



F(<Q)-Y,F{Qk 
k=l 



<MC<dCY^ ^dQk =--Y,h{Qk). 

LP{I) k=l k=l 



The function b : Fa^CZ'^) M.,Q i-^ dC^dQ satisfies the three conditions 
required for a boundary term. □ 
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Proof of Theorem [H Given the previous Lemma, we can apply the abstract 
ergodic theorem to our counting functions associated to finite subsets Q. In 
order to check the form of the error estimate, we note that the constant C 
in equation (fTTj) satisfies 

CP = C2 const(p, d) 

with C2 being the constant from Theorem [6] depending on the Kato-norm 
of Vc,-. Also, 

b{Q) = dC0Q 

so that the bound D on b in the general abstract ergodic theorem can be 
taken to be dC. It follows that ^^jg^ < Ad^C/M. 

Finally, the uniform lower estimate on the eigenvalues established in 
Lemma [5] gives a uniform upper estimate on the number of eigenvalues in 
(— 00, T], namely 

/ \ d./2 

Hence we can take 

K = C3{T + Ci)^ 

as the uniform bound on the function F in the abstract ergodic theorem. 
Here the constants Ci, C2, C3 depend only on d and the Kato norm of Vc,-- 

'□ 

5. Application to random operators 

In order to apply our results to random operators we quote the necessary 
random versions of the definition and main theorem from [GLVOSj : Let 
(r2,P) be a probability space such that acts ergodically on (J7,P). We 
denote the Z^'-action on Q hy x: uj t-^ u — x. A random ^-colouring is a 
map 

C : Q — > A with C{u - y)x~y = C{u})x 

for all X, y € Z"^. Note that for each fixed uj we obtain a (usual) ^-colouring. 
By the (usual) ergodic theorem for scalar functions the frequencies of pat- 
terns exist almost surely. Thus we can apply our abstract Banach space 
valued ergodic theorem: 

Theorem 13. Let A be a finite set, C be a random A-colouring and {X, || • ||) 
a Banach space. Let {Uj)j^fq be a van Hove sequence. For each fixed uj € 
let F^ : F'fin{1''^) — > X be a C{u!) -invariant, almost- additive bounded 
function. Assume that the family {F^)^^q is Z'^ -homogeneous, i.e. i*L+a;(Q+ 
x) = Fi^{Q) for all x G 7j'^,Q € J-fin{'^'^). Then, for almost every uj ^ Q the 
limits _ 

j^oo \Ui\ M^oo ^ \Cm\ 

exist in the topology of {X, \\ ■ ||) and are equal. In particular, F^ is almost 
surely independent of u. 
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Concretely, we apply this to colourings given by local models for the 
potentials like in Section [2l but now with a randomly chosen colouring so 
that all operators and counting functions additionally depend on the random 
variable w. For the formulation of the theorem we introduce a distribution 
function : M — ?> M defined by a trace per unit volume formula (sometimes 
called Pastur-Shubin formula) 



(12) 



N{X) := [ Tr[xH/oX(-oo,A](^c.)]dF(a;). 
Jn 



By applying the random version of the ergodic theorem we get the random 
version of Theorem [2) 

Theorem 14. Let C : ft — > (^^d A be a random colouring, ([/j)jgN a yo-n 
Hove sequence. Let / C M 6e a finite interval and p € [l,oo). Then for 
j ^ oo we have for almost all uj 



(13) 



1 



d\ 0. 



Remark 15. Of course there are similar error estimates as in Theorem [2j 



Proof. The existence of the limit in Equation (jl3p follows directly from our 
abstract theorem since we checked all requirements in the proof of Theorem[2] 
already. 

For the proof of the Shubin-Pastur formula ()12p we use a variation of the 
proof of Theorem 3 in |GLV07j : First notice that 



N{X) 



T7T I Tr[x^t/,X(-oo,A](^-)] dP(^) 



independently of Uj due to additivity and invariance. Now it suffices to 
show that 



(*) 



1 



■Tr 



-tHj 



for all t. We use the abbreviation 
of the trace to conclude 



rrUj 
llui 



0, j 



oo 



and the linearity 



Tr 



Tr [Xw^, (e" 



)] +Tr 



The second term actually vanishes. Thus it suffices to estimate the first 
term, which we do next. Let us note that for a compact operator K and a 
bounded operator B on the same Hilbert space the corresponding singular 
values obey the relation 

llniBK) < \\B\\fin{K). 

In particular this holds for K = T4ff as in Theorem [H] and B = X\Y^j • Now 
we can proceed exactly as in Lemma [T21 using that the addition of the 
boundary conditions by which ©if^^ differs from gives a boundary term 
in LP {I) for the corresponding SSF. Now the claim Q follows from the van 
Hove property. □ 
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